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MONOTONICITY AND CONDENSATION IN HOMOGENEOUS STOCHASTIC 

PARTICLE SYSTEMS 

THOMAS RAFFERTY, PAUL CHLEBOUN, AND STEFAN GROSSKINSKY 


Abstract. We study stochastic particle systems that conserve the particle density and 
exhibit a condensation transition due to particle interactions. We restrict our analysis 
to spatially homogeneous systems on finite lattices with stationary product measures, 
which includes previously studied zero-range or misanthrope processes. All known 
examples of such condensing processes are non-monotone, i.e. the dynamics do not 
preserve a partial ordering of the state space and the canonical measures (with a fixed 
number of particles) are not monotonically ordered. For our main result we prove that 
condensing homogeneous particle systems with finite critical density are necessarily 
non-monotone. On finite lattices condensation can occur even when the critical density 
is infinite, in this case we give an example of a condensing process that numerical 
evidence suggests is monotone, and give a partial proof of its monotonicity. 


1. Introduction 

We consider stochastic particle systems which are probabilistic models describing 
transport of a conserved quantity on discrete geometries or lattices. Many well known 
examples are introduced in [1], including zero-range processes and exclusion pro¬ 
cesses, which are both special cases of the more general family of misanthrope pro¬ 
cesses introduced in [2] . We focus on spatially homogeneous models with stationary 
product measures and without exclusion restriction, which can exhibit a condensation 
transition that has recently been studied intensively. 

A condensation transition occurs when the particle density exceeds a critical value 
and the system phase separates into a fluid phase and a condensate. The fluid phase is 
distributed according to the maximal invariant measure at the critical density, and the 
excess mass concentrates on a single lattice site, called the condensate. Most results 
on condensation so far focus on zero-range or more general misanthrope processes in 
thermodynamic limits, where the lattice size and the number of particles diverge simul¬ 
taneously. Initial results are contained in [3-5], and for summaries of recent progress 
in the probability and theoretical physics literature see e.g. [6-8]. Condensation has 
also been shown to occur for processes on finite lattices in the limit of infinite den¬ 
sity, where the tails of the single site marginals of the stationary product measures 
behave like a power law [9]. In general, condensation results from a heavy tail of 
the maximal invariant measure [10], and so far most studies focus on power law and 
stretched exponential tails [11,12]. As a first result, we generalize the work in [9] 
and provide a characterization of condensation on finite lattices in terms of the class of 
sub-exponential tails that has been well studied in the probabilistic literature [13-16]. 
This characterization holds for a particular definition of condensation given in Section 
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2.1, which was also used in [9]. Our main result is that all spatially homogeneous 
processes with stationary product measures that exhibit condensation on finite lattices 
with a finite critical density are necessarily non-monotone. 

Monotone (attractive) particle systems preserve the partial order on the state space 
in time, which enables the use of powerful coupling techniques to derive rigorous re¬ 
sults on large scale dynamic properties such as hydrodynamic limits (see [17] and 
references therein). These techniques have also been used to study the dynamics of 
condensation in attractive zero-range processes with spatially inhomogeneous rates 
[18-22], and more recently [23,24] . As we discuss in Appendix A, non-monotonicity in 
homogeneous systems with finite critical density can be related, on a heuristic level, to 
convexity properties of the canonical entropy. For condensing systems with zero-range 
dynamics, it has been shown that this is related to the presence of metastable states, 
resulting in the non-monotone behaviour of the canonical stationary current/diffusivity 
[25] . This corresponds to a first order correction of a hydrodynamic limit leading to an 
ill-posed equation with negative diffusivity in the case of reversible dynamics. Heuris- 
tically this is of course consistent with the concentration of mass in a small, vanishing 
volume fraction, but poses great technical difficulties to any rigorous proof of hydro- 
dynamic limits for such particle systems. First results in this direction only hold for 
sub-critical systems under restrictive conditions [26], and due to a lack of monotonic¬ 
ity there are no results for non-reversible dynamics. 

Condensing monotone particle systems would, therefore, provide interesting exam¬ 
ples of homogeneous systems for which coupling techniques could be used to derive 
stronger results on hydrodynamic limits. However, our result implies that this is not 
possible for condensing models with stationary product measures and a finite critical 
density on finite lattices. In the thermodynamic limit condensation has been defined 
through the equivalence of ensembles, which can be established in generality for a 
class of long-tailed distributions with a finite critical density [7, 10]. This class has also 
been studied before [14,15] and includes the class of sub-exponential distributions, 
for which our results apply also in the thermodynamic limit. A detailed discussion of 
their connections and the resulting differences between condensation on finite lattices 
and in the thermodynamic limit is given in Sections 4.1 and 4.2. We remark that for 
systems where the dynamics is directly defined on infinite lattices there are no rigorous 
results or characterizations of condensation to our knowledge, and we do not discuss 
this case here. 

For systems with infinite critical density condensation can still occur on finite lattices, 
and since non-monotonicity typically occurs above the critical density, such processes 
can also be monotone. When the tail of the stationary measure is a power law and 
decays faster than with the occupation number re, we prove that the process is 

still non-monotone. In Section 5 we present preliminary results for tails that decay 
slower than re“^/^, which strongly suggests that a monotone and condensing particle 
system exists (see [27] for further discussion). 

The paper is organised as follows. In Section 2 we introduce the background used to 
study condensation and monotonicity in particle systems, and state our main results. 
In Section 3 we prove our main theorem by induction over the size of the lattice, show¬ 
ing that the family of canonical stationary measures is necessarily not monotonically 
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ordered in the number of particles. In Section 4 we discuss the differences between con¬ 
densation on fixed lattices and in the thermodynamic limit, and prove equivalence of 
condensation on finite lattices with the tail of the maximal invariant product measure 
being sub-exponential. In Section 5 we review examples of homogeneous processes 
that have been shown to exhibit condensation, and present some explicit computations 
for misanthrope processes and processes with power law tails. 


2. Notation and results 

2.1. Condensing stochastic particle systems. We consider stochastic particle systems 
on fixed finite lattices A = {1,... ,L}, which are continuous-time Markov chains on 
the countable state space = N^. For a given configuration ij = {rj^ : x ^ A) E 
fli the local occupation, rj^ for x G A, is a priori unbounded. The jump rates from 
configuration 77 G Ol to ^ G are denoted 0 ( 77 , ^ 0, and the dynamics of the 

process is defined by the generator 

c(77,0(/(0 -/(77)) , ( 2 . 1 ) 

for all continuous functions / : —)■ M. We assume the process conserves the total 

number of particles 

Siiv) := ^ 77 ^ , 
xGA 

and conditioned on 5^ = N, the process is assumed to be irreducible, so that 5^ is 
the only conserved quantity. The process is therefore ergodic on the finite state space 
^L,N = {77 G Ol : 5 ^( 77 ) = W} for all fixed N ^ 0. On ^l,n the process has a unique 
stationary distribution and the family : At ^ O} is called the canonical 

ensemble. 

We focus on systems for which the stationary distributions are spatially homoge¬ 
neous, i.e. the marginal distributions vr^ ^^[ 772 ; G .] are identical for all x G A. This typ¬ 
ically results from translation invariant dynamics on translation invariant lattices with 
periodic boundary conditions, but the actual details of the dynamics are not needed 
for our results. For these systems we define condensation in terms of the maximum 
occupation number 

A/l( 77 ) := max 772 ; . ( 2 . 2 ) 

xeA 

Definition 2.1. A stochastic particle system with canonical measures ttl^n on ^l,n with 
L ^ 2 exhibits condensation if 

lim lim ttl n[Ml ^ N — K] = 1 . (2.3) 

if—>cx) V—>-oo ’ 

This condition implies in particular the existence of all limits involved. The inter¬ 
pretation of (2.3) is that due to the sub-exponential tails of the measure, in the limit 
W —)• 00 all but finitely many particles concentrate in a single lattice site and that 
the distribution of particles outside the maximum is non-degenerate. As we will see 
in Proposition 2.3 below, the latter is in fact given by the maximal invariant measure 
when the system exhibits product stationary measures. 
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There are of course other possible definitions of condensation which are less restric¬ 
tive or more appropriate in other situations. For inhomogeneous systems or thermody¬ 
namic limits (with N,L ^ oo) the condensed phase can be localized in particular sites 
or have a more complicated spatial structure (see e.g. [28] and for recent summaries 
[6-8] ). For our case of spatially homogeneous systems on finite lattices, A with fixed 
L, the above is the most convenient definition and has been used in previously studied 
examples [9]. A more detailed discussion is provided in Section 4. 


2.2. Monotonicity and product measures. We use the natural partial order on the 
state space given by rj ^ (ii and only if r]j; ^ Cx for all x G A. A function / : ^ M 

is said to be increasing if rj ^ ( implies /(r/) ^ /(C)- Two measures on Ql are 
stochastically ordered with /r ^ if for all increasing functions / : 17/, —> M we have 
//(/) ^ fJ-'if), where //(/) denotes the expectation of / with respect to p. 

A stochastic particle system on 17/, with generator £ and semi-group (S'(t) = : 

7 ^ 0) is called monotone (attractive) if it preserves stochastic order in time, i.e. 

fj. ^ fj.' ^ for all 7^0. 

Coupling techniques for monotone processes are an important tool to derive rigorous 
results on the large scale dynamics of such systems such as hydrodynamic limits. There 
are sufficient conditions on the jump rates (2.1) to ensure monotonicity for a large class 
of processes (see e.g. [17] for more details), however for our results we only need a 
simple consequence for the stationary measures of the process. 

Lemma 2.2. If the stochastic particle system as defined in Section 2.1 is monotone, then 
the canonical distributions tt/, tv tfre ordered in N, i.e. 

'^L,N ^ '^L,N+i for allN ^0 . (2.4) 

The proof is completely standard but short, so we include it for completeness. 

Proof Fix a monotone process as defined in Section 2.1. Consider two initial distribu¬ 
tions p and p', concentrating on 17 /,^tv and 17 /,^tv-i-i respectively, given by 

p[rj] = l{r]i = N) and = l(^i = + 1) , 

for p G 17/, TV and ^ G flL,N+i- Clearly p ^ p', and so by monotonicity of the process 
this implies pS{t) ^ p'S{t) for all 7^0. Furthermore, by ergodicity we have 

ttl,n = lim pS{t) ^ lim p'S{t) = TrL,N+i ■ 


All rigorous results on condensing particle systems so far have been achieved for 
processes which exhibit stationary product measures, for which the measures vt/t^tv 
take a simple factorized form. These can then be expressed in terms of un-normalized 
single-site weights w{n) > 0, n G N. Due to conservation of Sl such processes exhibit 
a whole family of stationary homogeneous product measures 

= n marginals f"'" ■ 

xGA 


(2.5) 
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The measures are defined whenever the normalization 

OO 

z{4>) := (2.6) 

n=0 

is finite. This is the case for all fugacity parameters (p G where = [0, pc) or 
[0,(^c], and 

pc '■= (limsup ^w{n)) ^ (2.7) 

n—>-cx) 

is the radius of convergence of (2.6). The family : p G £></>} is also called the 
grand-canonical ensemble and z{p) the (grand-canonical) partition function. Since the 
process is irreducible on ^Il.n for all G N we have w{n) > 0 for all n ^ 0. The 
canonical distribution can be written as 

= N] for all (/) G , (2.8) 

which is independent of the choice of p. Equivalently 

T^L,N['n] = II '^('^^) where Zl,n = JI w{r]^) (2.9) 

xgA V^^l,n x£A 

is the (canonical) partition function. Note that throughout the paper we characterize 
all measures by their mass functions since we work only on a countable state space 17 l 
and the measures 7rL,N concentrate on finite state spaces which is the framework 
we rely on in this paper. 


2.3. Results. Our results hold for systems with general stationary weights, w{n) > 0 
for each n G N, subject to the regularity assumption that 

lim w{n — l)/w{n) G (0, oo] (2.10) 

n—>-oo 

exists, which is then necessarily equal to pc- pc < oo then weights that satisfy (2.10) 
are sometimes called long-tailed [29], which is discussed in more detail in Section 4.2. 

For processes with such stationary product measures there is a simple equivalent 
characterization of condensation which we prove in Section 4.3. 


Proposition 2.3. Consider a stochastic particle system as defined in Section 2.1 with 
stationary product measures as defined in Section 2.2 satisfying regularity assumption 
(2.10). Then the process exhibits condensation according to Definition 2.1 if and only if 

Pc < oo, Dp, = [Q, Pc], and 


^Uvi+V2 = N] 

— -F- ITTi - 

N^oo up,^[r]i=N\ 


Z2. 


N 


lim , , , , 

N-^oo w{N)z[pc) 


G (0,oo) 


exists . 


( 2 . 11 ) 


In this case, the distribution of particles outside of the maximum converges weakly (equiv¬ 
alently in total variation) to the critical product measure i-e. for fixed ni,..., n^-i 
we have 


L-l 

TTL,N['ni = ni,...,pL-i=nL-i\ML = PL]^ Y\.i 2 p,fpi = ni] as N ^ oo . (2.12) 

i=l 
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Note that for (pc £ (0, oo) we may rescale the exponential part of the weights to get 
(/>c = 1 and we can further multiply with a constant, so that in the following we can 
assume without loss of generality that 


rc( 0 ) = 1 and pc = li™ w{n — l)/w{n) = 1 . 

n—>-oo 


The condition (2.11) can also be written as 


Z2,N 
lim — 
N^oo w{N) 


('w*w){N) 

lim --- 

N^oo W[N) 


G (0, oo) 


exists , 


(2.13) 


(2.14) 


where {w * w){N) = J2k=o w{k)w{N — k) is the convolution product. This is a standard 
characterization to define the class of sub-exponential distributions (see e.g. [30,31]). 
Sub-exponentiality implies that a large sum of two random variables is typically real¬ 
ized by one of the variables taking a large value (see Section 4.1 for more details), 
which is of course reminiscent for the concept of condensation. Implications and sim¬ 
pler necessary conditions on w{n) which imply (2.14) have been studied in detail, 
and we provide a short discussion in Section 4. As a consequence of Proposition 2.3, 
condensation is only a property of the tail of w, therefore if a process with stationary 
product measures (2.9) condenses in the sense of Definition 2.1 for some L ^ 2, it 
condenses for all L ^ 2 . 

Proposition 2.3 provides a generalization of previous results on condensation on 
finite lattices [9] and is used in the proof of our main result, which is the following. 


Theorem 2.4. Consider a spatially homogeneous stochastic particle system as defined 
in Section 2.1 which exhibits condensation in the sense of Definition 2.1, has stationary 
product measures that satisfy (2.10), and has finite critical density 

1 °° 

Pc ■= = -Tj-T nw{n)p'^ < oo , (2.15) 

n=0 

then the process is necessarily not monotone. 

The same is true if (2.15) is replaced by the assumption that^ w{n) ~ n~^ with b G 
(3/2,2]. 

2.4. Discussion. The class of distributions which fulfil (2.11) (called sub-expontial), 
and therefore exhibit condensation on finite lattices, is large (see e.g. [13, Table 3.7]), 
and includes in particular 

• power law tails w{n) ~ n~^ where b > 1, 

• log-normal distribution 

w{n) = - exp{-(log(n) - /i)^/( 2 (T^)} , (2.16) 

n 

where /r G M and cr > 0, which always has finite mean, 

• stretched exponential tails w{n) ~ exp{—C'n'*'} for 0 < 7 < 1 , C > 0 , 

• almost exponential tails w{n) ~ exp | — for /3 > 0 . 


^For functions /i : N ^ R we use the notation g{n) ~ h{n) if —>■ c G ( 0 , 00) as n —>• 00. 
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Density p = N / L 


Figure 1 . Non-monotone behaviour of the expected background den¬ 
sity E^f{N) ( 2 . 17 ) for lattice sizes L = 32 and L = 128; (A) (finite 
mean) power law tails with 6 = 5, (B) log-normal tails with /z = 0 and 
fj = l/-\/2, and (C) almost exponential tails with /3 = 1. The dotted 
black line shows the limit as L,N ^ oo and N/L ^ p, which is mono¬ 
tone and non-decreasing. 


For the last two examples all polynomial moments are finite. This covers all previously 
studied models on condensation in zero-range processes [5,9,12]. As we will discuss 
in Section 4.2 all these examples also exhibit condensation in the thermodynamic limit, 
whenever they have finite first moment. It can also be shown that the limit in (2.14) is 
necessarily equal to 2z{(f)c) and that in fact ^ Lz{(j)c)^~^ for any fixed L > 2 (see 

[16] and Proposition 3.2). 

Since we consider a fixed lattice A, pc < oo is not a necessary condition for conden¬ 
sation as opposed to systems in the thermodynamic limit. Even if the distribution of 
particles outside the maximum has infinite mean, condensation in the sense of Defini¬ 
tion 2.1 can occur. However, if z{(f)c) = oo (e.g. for power law tails with b ^ 1), the 
distribution of particles outside the maximum cannot be normalized, condition (2.11) 
fails, and there is no condensation in the sense of our definition. 
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We will prove non-monotonicity in the next section by showing that expectations 
for a particular monotone decreasing observable / : ^ M under 7tl,n are not 

decreasing in N. The chosen function is related to (but not equal to) the number 
of particles outside the maximum (condensate), which has been shown previously to 
exhibit non-monotone behaviour for a class of condensing zero-range processes in the 
thermodynamic limit [12,25]. When the number of particles N > pcL just exceeds 
the critical value, typical configurations still appear homogeneous with a maximum 
occupation number of^ on{N). Only when the number of particles is increased further 
the system switches to a condensed state with a maximum that contains a non-zero 
fraction of all particles. We present numerical evidence of this non-monotone switching 
behaviour for the background density 

:= - Ml) (2.17) 


in Figure 1. This is a finite size effect which disappears in the limit L —^ oo, and for 
specific models it has been shown to be related to the existence of super-critical homo¬ 
geneous metastable states [12,25,32]. For large L, the switching to condensed states 
occurs abruptly over a relatively small range of values for N. Since the ttl^n are con¬ 
ditional product measures the correlations in the system are very weak, which causes 
metastable hysteresis effects and non-monotonicity of the canonical ensemble around 
the critical point. Metastable hysteresis has been established in [12,25] for zero-range 
processes. Our result implies that this behaviour is generic for all condensing systems 
with finite critical density. We also give a heuristic discussion of the connection to 
convexity properties of the entropy of the system in Appendix A. 

There are several examples of homogeneous, condensing, monotone particle systems 
with finite critical density which have been studied on a heuristic level and which we 
summarize in Section 5. Their stationary measures are not of product form and no 
explicit formulas are known, so these systems are therefore hard to analyse rigorously. 
For systems with non-product stationary measures, upward fluctuations in the density 
which are homogeneously distributed may be suppressed strongly enough, so that the 
metastable states do not exist. Such models may then also be monotone, and examples 
are given in Section 5.3. 

We excluded the case </>£ = 0 in the presentation in Section 2.3 for notational conve¬ 
nience, but it is easy to see that our results also hold in this case. With the convention 
0° = 1 we have ^;(0) = w{0) = 1 and pc = 0, and then existence of the limit Z 2 ^n/w{N) 
is equivalent to 


'^2,n[M2 = A^] = 2 


w{N)w{0) 

Z2,N 


2w{0) 

2z(0) 


= 1 as N ^ oo, 


i.e. condensation of all N particles on a single site. This can easily be extended to all 
L ^ 2 with Proposition 3.2. Considering only events with all N particles on one site, or 
— 1 particles on one site and 1 particle elsewhere, we have convergence from above 


Zl,n 

w{N) 


Lw{0)^-^ > L{L - > 0 . 

w[N) 


^For functions /, g : N —>• R we use the notation g — On{f) if ^ 0 as n —>■ oo. 
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This implies the non-monotonicity of 7rL,N as discussed in Section 3. 


3. Proof of Theorem 2.4 


We assume that the process exhibits condensation in the sense of Definition 2.1 and 
has stationary product measures, so the canonical measures ttl^n are of the form (2.9). 
Furthermore, we assume the weights satisfy the regularity assumption, and without 
loss of generality 4>c = 1, see (2.13). We show that the family of canonical measures 
is not stochastically ordered in N, which implies non-monotonicity of the process by 
Lemma 2.2. To achieve this, we use the test function 

/(r/) := l(r/i = r ]2 = ... = t/l_i = 0) , (3.1) 

which indicates the event where all particles concentrate in the maximum at site L. 


Lemma 3.1. The function / : —)• M defined in (3.1) is monotonically decreasing, 

which implies that 


Zl,N , 


for all N ^ 0 , 


w{N) w{N + l) 

whenever the family of canonical measures 'Kl,n is stochastically ordered in N. 


(3.2) 


Proof Fix configurations pX G Dl such that jy ^ If /(jy) = 0 then p has at least one 
particle outside of site L, therefore so does ( which implies f{() = 0.lif{p) = l then 
necessarily f{p) ^ f{Q since f{() G {0,1}. Therefore / is a decreasing function. Using 
(2.9) and the convention (2.13), we find that the canonical expectation of the function 
(3.1) is given by 

w{0)^~^w{N) w{N) 

^ - • 

So if the canonical measures are monotone in N, monotonicity of / implies (3.2). □ 


By Proposition 2.3 we know that for condensing systems the ratio Z 2 ^n/w{N) con¬ 
verges. This implies the sequence /w{N) in Lemma 3.1 converges, as summarised 
in the following proposition using our notation (see [16, Theorem 1 and Lemma 5]). 


Proposition 3.2. Consider conditional product measures (2.9) with weights w{n) > 9 for 
all n G N, which satisfy 


w(n—l) I 

/ ^ Qr as n 

w{n) 

z{4>c) < OO, 

Z2,N 


OO, the regularity assumption (2.13), 


' TJiN) ^ C as N ^ OO. 
Then C = 2z{(j)c) and furthermore. 


Zl,n 

w{N) 


Lz{(j)c)^ ^ 


as N ^ OO for all L ^ 2 . 


(3.3) 


Note that the limit in (3.3) states that the probability of observing a large total 
number of particles under the critical product measure is asymptotically equivalent 
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(A) 


(B) 



Density p = N / 2 



Figure 2. Non-monotone behaviour of Hl{N) (3.4), which is the ex¬ 
pected value of the observable (3.1) rescaled by its limit; (A) power law 
tails w{n) ~ n~^ for L = 2 with 6 = 3, 1.75 and 1.35, where the latter is 
conjectured to be monotone (see Section 5); (B) log-normal tails (2.16) 
with /I = 0 and a = l/-\/2 for L = 2, L = 32 and L = 128. 


to the probability of observing a large number of particles on any one of the L sites, 
precisely 

< [Sl{v) = , 

1™ -7^-F-= 1 • 

Lu^^[r]i = N\ 

This is further equivalent to the canonical probability of the maximum containing the 
total mass converges to the critical probability that L — 1 sites are empty, i.e. til ,n[Ml = 
= 0 ]- 

To complete the proof we show that a subsequence of /w{N) converges from 
above, which contradicts the assumption of monotonicity by Lemma 3.1. We present a 
numerical illustration for the monotonicity properties of the function 


Hl{N) := 


1 

Lz(l)^-i w{N) 


(3.4) 


in Figure 2, which is normalized such that Hl{N) —)• 1 as —)■ oo. The proof of the 
following lemma represents the most significant part of the proof of Theorem 2.4 and 
is given in Section 3.1 for the case of finite mean and in Section 3.2 for the power law 
case. 


Lemma 3.3. Under the conditions of Theorem 2.4, and assuming without loss of general¬ 
ity fc = h for each L ^ 2 there exists a C > 0 and a sequence Nm G N with Nm oo as 
m —)• oo such that 


min 

71 ^ Nm 


ZL,n 

w{Nm) 



^ C/Nm , 


Therefore, we know that there exists some A^* G N such that Zl,n* /w{N*) > 
Zl,n*+i/w{N* + 1), which by Lemma 3.1 implies that the canonical measures are 
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not stochastically ordered in N, and thus the process cannot be monotone by Lemma 
2.2. This completes the proof of Theorem 2.4. 


3.1. Proof of Lemma 3.3: The finite mean case. In order to prove that con¬ 

verges from above for some non-decreasing sequence Nm, we first specify a sequence 
on which we can bound the ratio below. 

w{Nm) 


Claim 3.4. For weights {w{n) : n G N} with finite and non-zero first moment, i.e. 0 < 
Pc < oo, there exists a sequence Nm G N with Nm — )• oo as m — )• oo such that for all 
k G {0,... , Nm — 1} 


w{Nm - k) ^ . k 
w{Nm) '' Nr, 


(3.5) 


Proof For each m G N, define Nm as follows 


Nm = max{n ^ m : nw{n) = min j w{j)} . 

j 

By definition Nm is a non-decreasing sequence. Assume for contradiction that Nm is 
bounded above, then for all j G N we would have jw{j) ^ j*w{j*) > 0 for some j* G N, 
and therefore nw{n) —)• oo contradicting the assumption of finite first moment. For 
A: G {0,... Nm — 1} we have 


{Nm - k)w{Nm -k)^ NmW{Nm) 


and thus 


w{Nm - k) ^ Nm ^ Jf 
w{Nm) " (Nm-k)" ^ Nr 


□ 


Claim 3.5. For weights {w{n) : n G N} with finite first moment, there exists a subsequence 
{Ni : f G N} o/ the sequence defined in Claim 3.4 such that > 22(1) for all 

n G {0,..., Ni} and i sufficiently large. 

y 

Proof By neglecting at most a single term in the sum defining Z 2 ,n, the ratio can 
be bounded below as follows, 


Z2,N 

w{N) 


N 

'^w{n) 

n=0 


w{N — n) 
w{N) 


^ 2 


L7V/2J-1 

E 


n=0 


w{n) 


w{N — n) 
w{N) 


We define 


Km. := max <! k* ^ Nm '■ = min Z- 


0 s£ fc ^ V„ 


J2,fc 


to be the largest index where the ratio 


Z2,k 

w{Nm) 


is minimized. In particular 


m ^ 0 and n G {0,..., Nm}- 

w{Nm) w{Nm) 


(3.6) 


(3.7) 


By definition Km ^ Nm, and so r := Yrmswpm^oo^m/Km ^ 1- There exists a subse¬ 
quence {mi : i ft 0) such that Km^/Nm^ r, with a slight abuse of notation we denote 
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the subsequences and Kmg simply by and Kf. Suppose r < 1, by Claim 3.4 we 
have 


^ Z2,Ki 
w{Ni) '' w{Ki) 


2 - 


Ni 


2z(l)(2 -r) > 2z(l), 


which together with (3.7) contradicts Proposition 3.2, therefore Ki/Ni —)• 1 and K^/Ni ^ 1 
for all f. 

Applying Claim 3.4 we then have 


Z2,K^ 

w{Ne) 




^ 2 w{k) 


k=0 

Subtracting 2z(l) we get 
Z2,Ki 


w{Ki - k) 
w{N() 


{_Ki/2\-l 


> Y1 kw{k)+2(2 




[Ke/2\-l 


k=0 


NeJ 




k=0 


w{Ni) 


-2z{l) 

„ IK,/2\-1 oo , IK,/2\-1 

^ kw{k)-2 ^ y;(A:)+2[l- ^ 


k=0 fc=[A'f/2j 

Neglecting the final term in (3.8) we have 

^ LA'^/2J-1 

Z2,Ki 


Ni 


w{k) . (3.8) 


k=0 


Np 


w{Ni) 


2z(l) j >2 kw{k) — 2Np E 


fc =0 


k=[Kp/2\ 


\_Kf,/2\-l 


>2 Y^ kw{k) -^ Y^ kw{k) ^ 2pcz{l) > 0. 


k=0 


K, 


k=lKe/2\ 

as £ —)• OO, where pc is the critical density defined in (2.15). Together with (3.7) this 
completes the proof of Claim 3.5. □ 

To complete the proof of Lemma 3.3 we proceed by induction on the system size, L. 
We make the following inductive hypothesis; 

(H) there exists a sequence {Nm ; m G N} such that > Lz{l)^~^ for all 

n G {0,..., Nm} and m sufficiently large. 

The case L = 2 is given by Claim 3.5. Analogously to the proof of Claim 3.5 we define 


Km := max \k* Nm'- Zih* = min Zik t- 


(3.9) 


By the same argument as in the proof of Claim 3.5 there exists a subsequence {mi : 
£ ^ 0) such that Kmp/Nmp —^ 1, again we denote the respective subsequences by Ki 
and Np. For £ sufficiently large, we have 




Kp 


ZL,Ki-k 


VKp/2\ 


w{Ni) 

^ fc=0 ^ ’ k=0 


7 VKp/2\-l 

E p,^ZLKt-k ry w{Ki — k) 


LA'^/2J lKp/2\-l 

> L z{l)^ ^ w{k) + I 1 + 


k=0 


k=0 


k=0 

Ni-Ki + k 

Ni 


w{Ni) 
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where the final inequality follows from the inductive hypothesis (H) and Claim 3.4. 
Subtracting (L + we get 


w{Nf) 


OO 

{L-\-> — L ^ w{k) 

k=[Ki/2\+l 



[Ki/2\-l 


k=0 



lKe/2\-l 

+ ~ ^(1)^ • 


/c=0 


(3.10) 


Now, following the proof of Claim 3.5, multiply (3.10) by Ni and neglect the second 
term on the second line. Then the first term vanishes, since 

CXD 2>N ^ 

E < -!?7 E k w{k) —)■ 0 as f —>■ OO . 

fc=LA'^/2j+l ^ k=lKe/2\+l 


In terms of the normalized grand-canonical measure = k], so we 

have 

OO 

^ k ZL^k = z{1)^u^{Sl) = PcLz{l)^ G (0, oo) , (3.11) 

k=0 

where pc is the critical density as defined in (2.15). This implies that the first term in 
the second line of (3.10), after multiplication with Ni, converges to a strictly positive 
constant. Finally, the third line in (3.10) converges to zero after multiplying by since 
we have ZL,n = z{l)^, which implies 


1 , oo ^ oo 

0^iV,( ^ ZL^k-z{l)^] =-N, E kZL,k^0 

^ fc=o ^ fc=[A:^/2j fc=[ir^/2j 

as f oo, by (3.11). Using the definition of Ki in (3.9), this implies that there exists a 
constant (7 > 0 such that for all i large enough 


min 

n ^ N£ 


f ZL-\-l,n 

V wi^i) 


(L + 1)^(1)^^ 


^ C/Ne , 


so (H) holds for L + 1, completing the induction. This concludes the proof of Lemma 
3.3 for the case where the critical measure has finite mean. 


3.2. Proof of Lemma 3.3: The infinite mean power law case. We consider stationary 

weights of the form w{n) = n' ^h{n) with ri;(0) = 1, h{n) ^ c G (0,oo), and b G (1,2). 
We prove non-monotonicity of Z^^^/w{N) for b G (3/2,2) and h{n) = 1 for all n G N 
via an exact computation. The case 6 = 2 can be done completely analogously but 
involves different expressions with logarithms in the resulting limits, and is presented 
in Appendix B. The proof remains valid for general converging h{n) with only minor 
differences, which we explain in a remark at the end of this section. Convergence of 
Z 2 ,n/'w{N) —^ 2z(l) from above or below for the exact power law depends on the 
parameter 6 G (1, 2), as summarized in the next result. 
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Lemma 3.6. For stationary weights of the form w{n) = n ^ and w(0) = 1 with b G (1,2) 

N’’-^(^^-2zil)]^F2ib) asN^oo, (3.12) 

\w[N) f 

where 


OO .. 2—1 

^46)=2^-n(i+6) 


2 ^— 1—2 

- 2 - 


1=1 j=0 


1 — b + i 


b-1 



i/ 6 G (|, 2 ) 
if be ( 1 , |) 


For L > 2we have 

= FL{b) := z{l)FL-i{b) + {L- l)z{l)^-^F^ib) , 

(3.13) 


which has the same sign as ^ 2 ( 6 ). 

This result implies that whenever w{n) = n~^ for n ^ 1 and w{0) = 1 with b G 
(3/2,2) Lemma 3.3 holds with C = ^ 2 ( 6 ). This completes the proof of Lemma 3.3 in 
the case h{n) = 1 . 


Proof of Lemma 3.6. To prove this result we make use of the full Taylor series of (1 — 
x)“^ at X = 0 and integral approximations to compute the asymptotic behaviour of 
summations. To simplify notation we assume that N is even. For odd N there is no 
term with multiplicity one and there exists an obvious modification. First note that 
w{n) fulfils the regularity assumption (2.13) and /w{N) —)■ Lz{l)^~^ as N ^ 00 
for all L ^ 2 [9], so by Proposition 2.3 a process with stationary measures will 
exhibit condensation. For L = 2 we subtract 2^(1) from Z 2 ^n/'w{N) to get 


Z2,N 

w{N) 


2 z(l) 


N/2 


^ , ,w(N-n) ^ , w(N/2)w(N/2) 


71=0 

N/2 




n 


-b 


71=0 


1 - — 
N 


71=0 
—b 00 


w{N) 
2 ^n"^ - . 


n=0 


Substituting the Taylor expansion of (1 — x) ^ we find 


(3.14) 


Z2,N 

w{N) 


2 z(l) 


n /2 00 ^ i 1—1 00 

71=0 1=0 ■ ^ ^ j=0 


OO i—1 


N/2 


2E/n(2+'')4E"^“ 


22=1 

OO 


Z! 

2=1 j=0 


n=l 


2 . (3.15) 

n=Ar/2+l 


In the last line the i = 0 term was combined with the second term, and we adopt the 
usual convention that empty products are equal to one. Both summations in n are over 
continuous and monotone functions g : R ^ {0, 00), therefore we can use the usual 
integral approximation for decreasing (increasing) functions 



d 

g{x) dx din) ^ ^ )gic) + 


n=c 


g{x) dx 


(3.16) 
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for all c G N and d G N U {oo}. Multiplying by ^ we find the limit as M —)• oo of 
(3.15) to be 


CXD ^ 2—1 

=2E^n{i+ 

i=l j=0 

It is shown in Appendix C that this is positive (and finite) in the region b G (3/2,2) and 
negative (and finite) in the region b G (1,3/2), completing the proof of Lemma 3.6 for 
L = 2. The result holds for general system size, L ^ 2, and is proved by induction. The 
inductive h 5 ^othesis states 


b—l—i 


'fb-l 




- 2 - 


-b+i 6-1 


(3.17) 


lim N^' 

N^OO 


-1 f Zl,n 


-Lz{l)^-^ ) =Fl(6) = z{l)FL-iib) + {L-l)zil)^-'^F2{b) . (3.18) 


L-2 


w{N) 

Similar to the case L = 2 we write 


N' 


i-lf Zl+1, 


N 


\ w{N) 


= N‘ 


6-1 


N/2 




n=0 


w{N — n) 
' w{N) 


- z 




N/2-1 


(lY +iV‘-l w(n 


n=0 


ZL,N-n 

w{N) 


-LzilY 




(3.19) 


We first establish the limit of the function Ql,n in equation (3.19). The inductive 
hypothesis (3.18) can be written as 


ZL,n _ FL{b) + On(l) 


wyn) n 

which implies Ql,n can be written as 

^N/2-l 


6-1 


+ Lz(V) 


L-l 


Ql,n = N' 


6-1 


^ / .w{N — n) Zi^N-n J 

.../An FJTt - 


72=0 
^ N/2-1 


w{N) w{N — n) 


YY, w{n) 


w{N — n) 


72=0 




w{N) [ {N - nf-^ 


(3.20) 



Rearranging terms and noting that ^ have 

/AT — \ 

©L.V =(-F’l( 6)+ 07v(l)) Y 

n=0 ^ ^ 

.AI/2-1 \ 

+ -(!))■ 

After Taylor expanding (1 — appearing in the first line above, it is easy to see that 

the limit of the first line is given by Fl(6)z( 1) as M —)■ oo. Using the L = 2 result to 
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calculate the limit of the second line we find 

Ql,n FL{b)z{l) + ^ a,s —)• oo . (3.21) 

To identify the limit of Hx, tv in (3.19), we again make use of the Taylor expansion of 
(1 — x)“^ similarly to the two site case and we write 

, N/2 oo ^ i-1 / \ * \ 

Hl,TV = ^ ^ n (■?■ + • 

^ n=0 1=0 ■ j=0 ^ ^ ^ 

Changing the order of summations, separating the i = 0 term and using ^L,n = 

z{l)^ we have 

y oo i—l N/2 oo \ 

= x: Jill O'+ '')«! • ».22) 

^ i=l j=0 n=l n=Af/2+l ^ 


For alH ^ 1 and b € (1,2) we have N’^ ^ > 0 as —)■ oo, which implies that for any 
fixed A^i e N we have n*Zx,n —^ 0. Therefore, the following limits are 

equal 


lim El,n = 

N—^oo N—^oo 


oo - 2—1 


N/2 


' i=l j=0 n=N\ n=N/2+l 

Using the inductive hypothesis (3.20) we have limTv-s-oo '^l,n is given by 

y OO i-1 n/2 

V>-TF,(6) + „„(!)) ( E IT n U + '>) E “• 


l\ 

l=l j={) 
oo ^ 2—1 


.r- / ^ 

at* n' 

n=Ni 


6-1 


E 


w{n) 


n=N/2+l 


n 


b-l 


N/2 


/ ^ V — j. - / 

+(E X n (•?■+^) ]^ E - Y. « 

n=Ni n=N/2+l 


i=l j=0 


Now applying the L = 2 result it is possible to show that 

Lz{l)^-^F2{h) 

^L,N —^ -X- , 


(3.23) 


(3.24) 




- {L + l)z(l)'* 


2 :(l)Fx( 6 ) + Lz(l)'^ ^F 2 {b) as A^ —)• oo , 


where the limit of the first line of (3.23) was 0 by the additional factor 1/n^ ^ appearing 
in the summations. Combining (3.21) and (3.24) we have 

jb-i ( Zl+i,n 

V wiN) 

concluding the induction so the result holds for all L ^ 2. From the recursion (3.13) it 
is obvious that Fiib) will have the same sign as ^ 2 ( 6 ), completing the proof of Lemma 
3.6. □ 


A slightly modified version of Lemma 3.6 also holds if the stationary weights are of 
the form w{n) = n~^h{n) where lim„_!.oo h(n) = c G (0, 00 ). The limit in (3.12) only 
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depends on the tail behaviour of the weights and is now given by cF 2 {b). Briefly, this 
can be seen as follows, (3.14) becomes 


NI2 

2 ^ n~^h{n) 
n=0 


h{N — n) 

h(N) 



2 V „-/.(-)+ 

^ h{N) 

n=0 ^ ’ 


Taylor expanding (1—x) ^ and rearranging to And terms of the form ^ * Y^=\ h{n)n 
and using the same argument to calculate the limit of '^l,n we have 


N/2 N/2 

lim y h{n)n-^+^ = lim V cn-^+^ < oo 

N^oo N^oo 

n=l ri=Ni 

for all i ^ 1 and any Ni e N, and the result follows. Similar modifications are required 
in the inductive step and the new limit in (3.13) is given by c^~^FL{b) for all L ^ 2. 
This does not change the sign of the limit in (3.13) and therefore Lemma 3.3 still holds. 


4. Characterization of condensation 


Condensation arises in spatially homogeneous systems with stationary product mea¬ 
sures due to the sub-exponential tail of the stationary weights w, which has been stud¬ 
ied extensively in previous work. In this section we review relevant results on heavy¬ 
tailed distributions and discuss the links between condensation on finite fixed lattices 
and in the thermodynamic limit before we give the proof of Proposition 2.3 in Section 
4.3. 

4.1. Sub-exponential distributions. Sub-exponential distributions are a special class 
of heavy-tailed distributions, the following characterization was introduced in [33] 
with applications to branching random walks, and has been studied systematically in 
later work (see e.g. [14-16,30]), for a review see for example [13] or [31]. 

A non-negative random variable X with distribution function F{x) = P[X ^ x] is 
called heavy-tailed if F{0+) = 0, F{x) < 1 for all x > 0, and 

gAx(i _ OO as X —>■ OO for all A > 0 . (4.1) 


It is called sub-exponential if F{0+) = 0, F{x) < 1 for all x > 0, and 

1 _ F*^(x) 

- -^2 as X ^ oo . (4.2) 

1 — F[x) 

Here F*^(x) = P[2fi + X 2 ^ x] denotes the convolution product, the distribution 
function of the sum of two independent copies Xi and X 2 . It has been shown [33,34] 
that (4.2) is equivalent to either of the following conditions. 


1 _ F*^(r) 

lim --V = L for all L^2, or 

x^oo 1 - F{x) 


lim 




x-i-oo p[max{Xj : i G {1,..., L}} > x] 


= 1 for all L ^ 2 . 


(4.3) 

(4.4) 


The second characterization shows that a large sum of independent sub-exponential 
random variables Xi is typically realized by one of them taking a large value, which 
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is of course reminiscent of the condensation phenomenon. It was further shown in 
[13,33] that sub-exponential distributions also have the following properties, 

lim ^ 7 ^ =1 Vy e M, (4.5) 

x^oo 1 — [X) 

e“(iF(x) = oo Ve > 0 (no exponential moments), (4.6) 

F{x)e'^^ —)■ oo Ve > 0 (slower than exponential decay) . (4.7) 



Most results in the literature are formulated in terms of distribution functions and tails 
and apply to discrete as well as continuous random variables. [16] provides a valuable 
connection to discrete random variables in terms of their mass functions ui(n), n e N. 
Assume the following properties for a sequence {w{n) > 0 : n € N}, 

1 as n OO, 

(b) z{l) := ^ (Ojoo) (normalizability), 

(c) limAT^oo = C G (0, oo) exists. 

Then [16, Theorem 1] asserts that C = 2z{l) and w{n)/z{l) is the mass function of a 
discrete, sub-exponential distribution. The implication 

Lz{l)^-^ as iV ^ oo for L > 2 
w{N) ^ ’ 

is given in [16, Lemma 5]. Sufficient (but not necessary) conditions for assumption (c) 
to hold are given in [16, Remark 1]. 

Provided z{l) < oo, then (c) holds if either of the following conditions are met: 

© win—k) ^ 

SUPi ^ fc ^ n/2 ^ ^ 

for some constant AT > 0, or 
(ii) 'u;(n) = e“”'V(”') 

where ?/)(x) is a smooth function on M with 'ijj{x) \ 0 and x‘^\'ilj'{x)\ ^ oo as x ^ oo, 
and < oo. 

Case (i) includes distributions with power law tails, w{n) ~ n~^ with 6 > 1. The 
stretched exponential with 'ijj{x) = x'^~^, 7 G (0,1), and the almost exponential with 
V’(x) = (log(x))“'^, /3 > 0, are covered by case (ii). The class of sub-exponential distri¬ 
butions includes many more known examples than the list given in Section 2.4 (see e.g. 
[13, Table 3.7]). Analogous to the characterisation of sub-exponential distributions, 
given by (4.4), for discrete distributions the existence of the limit {w * w){N)/w{N) is 
equivalent to the existence of the following condition 


P[Xi + X2 = N] 
P[max{Xi,X2} = N] 


—>■ 1 as —)• 00 . 


(4.8) 


This holds, since we have the following equality of ratios 


p[a:i + a :2 = iv] _ Z2,n _ {w*w){n) 

P[max{Xi,X 2 } = iV] “ 2w{N)Y:n=ow{n) ~ 2w;(iV) ’ 

Specific properties of power law tails w{n) are used in [9] to show condensation 
for finite systems in the sense of Definition 2.1. In Proposition 2.3, proved in Section 
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4.3, we extend this result to stationary product measures with general sub-exponential 
tails. In this context, condensation is basically characterized by the property (4.4) 
which assures emergence of a large maximum when the sum of independent variables 
is conditioned on a large sum. As summarized in the introduction, condensation in 
stochastic particle systems has mostly been studied in the thermodynamic limit with 
particle density p ^ 0, where L,N ^ oo such that N/L ^ p. In that case conditions 
on the sum of L independent random variables become large deviation events, which 
have been studied in detail in [35,36]. 


4.2. Connection with the thermodynamic limit. In the thermodynamic limit, a def¬ 
inition of condensation is more delicate and the approach presented in [7, 10] follows 
the classical paradigm for phase transitions in statistical mechanics via the equivalence 
of ensembles (see e.g. [48] for more details). A system with stationary product mea¬ 
sures (2.5) exhibits condensation if the critical density (2.15) is finite, i.e. pc < oo 
and the canonical measures 7rL,N are equivalent to the critical product measure in 
the limit L,77 —> oo such that N/L —^ p for all super-critical densities p ^ pc- The 
interpretation is again that the bulk of the system (any finite set of sites) is distributed 
as the critical product measure in the limit. It has been shown in [10] (see also [7] 
for a more complete presentation) that the regularity condition (2.10) and pc < oo 
imply the equivalence of ensembles, which has therefore been used as a definition of 
condensation in [7, Definition 2.1]. Therefore, any process that condenses for fixed A 
with Pc < oo in the sense of Definition 2.1 also condense in the thermodynamic limit. 
This includes all previously studied examples [5,9, 12], however there exists distribu¬ 
tions that satisfy (2.10) with pc < oo but do not satisfy the conditions of Proposition 
2.3 and do not condense for fixed A. This is illustrated by an example given below. It 
is also discussed in [7, Section 3.2] that assumption (2.10) is not necessary to show 
equivalence, but weaker conditions are of a special, less general nature and are not 
discussed here. Note also that equivalence of ensembles does not imply that the con¬ 
densate concentrates on a single lattice site, the latter has been shown so far only for 
stretched exponential and power-law tails with pc < oo in [11,12]. Both definitions 
involve only the sequence 7rL,N of canonical measures and not on the dynamics of the 
underlying process. Since the canonical measures (2.9) are fully characterised by the 
weights w{n) condensation can be viewed as a feature of the tails of the weights w{n). 

The condensation phenomena can also be studied for continuous random variables 
on the local state space [0,oo), see for example [36]. The following continuous ex¬ 
ample, taken from [15] is shown to satisfy (2.10) but is not sub-exponential. We 
show that the distribution has a finite mean and therefore exhibits condensation in 
the thermodynamic limit but not on a finite lattice in the sense of Definition 2.1. For 
a real-valued random variable X with distribution function F{x) = ^ x], assume 

F'{x) = g'{x)e~^^^\ Let (xn)„gH to be an increasing sequence with xq = 0 and g{x) 
be a continuous and piecewise linear function such that g{0) = 0 and g'{x) = 1/n for 
x e {xn-i,Xn)- The sequence (xn)ngN is defined iteratively as follows 
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Xn — Xn-1 = 

gixn) - g{xn-i) = , 


and g{x) — g{xn-i) = ^ ^ for x G [xn-i,Xn)- The mean is finite since 



where the final step uses the relation g{xn) — g{xn-i) = ^ 2(1 + g{xn-i)) to 

bound the series from above. 

For all distributions satisfying ( 2 . 10 ) which are not sub-exponential Zl,n/w{N) does 
not have a limit in (0, oo) as oo and with Proposition 2.3 there is no condensation 
on finite lattices according to Definition 2 . 1 . For a discretized version of the example 
given above with weights w{k) = exp{—g{k)} we have Z 2 ^n/w{N) —> oo for ~ as 
n —)• oo [15]. For this example, following the proof of Proposition 2.3 this implies that 
'^2,N[gi A r/2 ^ AT] —> 0 for A^ ~ Xn as A^ —> oo and all K ^ 0 . Therefore, the L = 2 bulk 
occupation number gi A t/2 diverges in distribution as A^ —> oo by receiving a diverging 
excess mass from the condensate due to the light tail of w{n). It can be shown that 
these weights satisfy (2.10) and therefore exhibit condensation in the thermodynamic 
limit where the excess mass can be distributed on a diverging number of sites. 

For a process that exhibits condensation in the thermodynamic limit with a sub¬ 
exponential critical product measure. Proposition 2.3 implies that condensation occurs 
also on finite lattices with pc < oo. Theorem 2.4 then implies that this process is neces¬ 
sarily non-monotone for all fixed system sizes L. However, monotonicity for condensing 
processes with long-tailed but not sub-exponential stationary measures remains open. 

4.3. Proof of Proposition 2.3. Let us first assume that the process exhibits condensa¬ 
tion according to Definition 2.1 and has canonical distributions of the form (2.9) where 
the weights fulfil (2.10), i.e. w{n — l)/w{n) cpc ^ (0, oo] as n —)• oo. In this part of 
the proof we establish that; 

(1) (t)c < oo, 

(2) has a limit as A^ —)■ oo, 

(3) z{(f)c) < oo, which also implies as A^ —> oo, and 

(4) convergence of —s- Lz{(l)c)^~^ for some L ^ 2 implies convergence for 

L = 2 and therefore (2.11) holds. 

Step (1), show (f)c < oo. Assume first that w{n — l)/w{n) — oo as n —)■ oo. For all 
AT G N and N > K we have 



^ L— - max 

Zl^N 0 ^ n ^ K 


{ZL-i,n) max {w{N - n)) . 


0 ^ n ^ K 
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Let re* = argmaxo ^ n s£ — n)) ^ K. The partition function is trivially 

bounded below by the event that site 1 takes — re* — 1 particles and the second site 
takes the remaining re* + 1 particles, i.e. 

Zl,n ^ w{0)^~‘^w{n* + l)w{N — re* — 1) . 


Therefore 

^ N - K] 


L 


K + 1 w{N-n*) 


max 


w{n* + 1) w{N — re* — l)o^n<iY 


(Zl- 


l,n) 


as iV —^ oo, which implies condensation cannot occur in the sense of Definition 2.1 
contradicting the initial assumption, therefore (pc < oo. 

Step (2), prove Zl^n/w{N) converges as ^ oo. By Definition 2.1 the limit 


QK ■= lim ^ N - K] , 

N^OO 

exists and > 0 for K sufficiently large. For N > K we have 

K 


N - K]= L—^^ ^ ZL-l,n- 


w{N — re) 
w{N) 


(4.9) 


(4.10) 


Since w{N — n)/w{N) cp'^, K is fixed, and ax > 0, (4.10) implies the convergence 

of Zl,n/w{N) as A^ —)■ oo. 

Step (3), prove z{(pc) < oo. By (2.3) we have —)• 1 as iT —)• cx3, taking the limit as 

A^ —)> oo of (4.10) this implies 


K 


OO 


lim 

K^oo 


Y, ZL-l,n<Pc = Y ZL-l,n€ < OO . 

n=0 n=0 


(4.11) 


Since we also have ZL-i,n(pc = this implies z{(pc) < oo. Using ax 1, 

(4.10) then also implies Zl^n/w{N) — )■ Lz{<pc)^~^ as At —s- cc. 

Step (4). We have seen above that condensation implies (pc < oo, z{(pc) < oo, and 
Zl,n/w{N) — )• Lz{<pc)^~^ as N ^ oo, then [37, Theorem 2.10] implies 


lim ^ 

N^oo W{N) 


‘iz{(pc) , 


completing this part of the proof. 

Now, let us consider a stochastic particle system with canonical distributions of the 
form (2.9) which fulfil (2.13) and (2.14) with (pc = I and z{l) < oo. We keep the 
notation for (pc = '^ general in the following to clarify the argument. It is immediate 
from Proposition 3.2, and remembering that we set wpS) = 1, that 


^l,nWl = N] = Lw{N)/Zl,n ^ z{PcY^-^'> > 0 . 


Then we have for all fixed K and N > K 

K 




^m[Ml^N-K] = lY 


w{N - n)ZL-l,: 


n=0 


Zl,n 


K 

'Y, ZL-l,n 

n=0 


w{N — re) Lw{N) 
w{N) Zl,n 


K 


E 

n=0 


ZL-l,nPc 

Z{pc)^~^ 


^(pSVl + • • • + riL-i ^ K) 
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as —)• oo. Since is a non-degenerate probability distribution, this implies that 

+ ... + riL-i ^ K) ^ 1 as K ^ oo, which is (2.3). 

To compute the distribution outside the maximum we get for fixed ni,..., and 
large enough N 


'n'L,N[rji=ni, ■ ■. ,r]L-i=nL-i\ML=r]L\ = 
1 


L7rL,N[ML = Vl] 


w{ni) ■ --winL-i) 


w{ni) ■ ■ ■ w{nL-i)w{N — ni— ... —n^-i) 
'^l,n[M^l=Vl] Zl,n 
w{N — ni — ... — riL-i) Lw{N) 


w{N) 


Zl,n 




\i,-i 


(4.12) 


as —)• oo. Here we have used that spatial homogeneity of the measure and asymp¬ 
totic uniqueness of the maximum according to (2.3) imply = ijl] \/L. This 

completes the proof of Proposition 2.3. 


5. Examples of homogeneous condensing processes 

In this section we review several stochastic particle systems that exhibit condensa¬ 
tion. By Theorem 2.4, if these processes are homogeneous and monotone with a finite 
critical density they do not have stationary product measures. To prove monotonicity 
for the examples mentioned below it is sufficient to construct a basic coupling of the 
stochastic process which preserves the partial order and particles jump together with 
maximal rate. For a definition of a coupling see [38] and for the statement of Strassen’s 
theorem linking stochastic monotonicity and the coupling technique see [39] . The steps 
to construct a basic coupling are outlined in [17]. 


5.1. Misanthrope processes and generalizations. Condensation in homogeneous par¬ 
ticle systems has mostly been studied in the framework of misanthrope processes [2, 
17]. At most one particle is allowed to jump at a time and the rate that this occurs de¬ 
pends on the number of particles in the exit and entry sites. The misanthrope process 
is a stochastic particle system on the state space H/, = defined by the generator 

^ 'r{r]^,riy)p{x,y){f{ri^^y) - f{ri)) . (5.1) 

x,yGA 

Here = v — Sx + Sy denotes the configuration after a single particle has jumped 
from site x to site y, which occurs with rate r{r]x,riy). The purely spatial part of the 
jump rates, p{x, y) ^ 0, are transition probabilities of a random walk on A. Such 
models are usually studied in a translation invariant setting with periodic boundary 
conditions, typical choices are symmetric, totally asymmetric or fully connected jump 
rates with p(x, y) = l/2{5y^x+i+5y,x-i), p{x,y) = 6y^x+i, or p{x,y) = {l-6y^x)/{L-l), 
respectively. 

Misanthrope processes include many well-known examples of interacting particle 
systems, such as zero-range processes [1], the inclusion process [40,41], and the ex¬ 
plosive condensation model [42]. It is known [2,27] that misanthrope processes with 
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translation invariant dynamics p{x, y) = q{x — y) exhibit stationary product measures 
if and only if the rates fulfil 


r(n, m) 

r(m + 1, n — 1) 

and, in addition, either 


r(re, 0)r(l, m) 
r{m + 1,0)r(l, n — 1) 


for all n ^ l,m ^ 0 , 


{ q{z) = q{—z) for all z £ A or, 
r(n, m) — r{m, n) = r{n, 0) — r(m, 0) for all n, m ^ 0 . 

The corresponding stationary weights satisfy 


w{k + l) w{l) r{l,k) -Pjr{l,k-1) 

w(k) ui(0) r(A: + 1,0) r(k,0) 

k=l 


(5.2) 


(5.3) 


(5.4) 


Misanthrope processes are monotone (attractive) [2] if and only if the jump rates satisfy 


r{n,m) ^ r{n + l,m) i.e. non-decreasing in n, 

r{n,m) ^ r(n,m-|- 1) i.e. non-increasing in m . (5.5) 


In Theorem 2.4 we have proved that processes that exhibit stationary product mea¬ 
sures and condensation with finite mean or power law tails, w{n) ~ n~^, with b G 
(3/2,2] are necessarily not monotone. For power law tails with b G (1,3/2] conver¬ 
gence of Zl,n/w{N) is from below and our method does not disprove monotonicity of 
the measures 7tl,n or monotonicity of the underlying process. Using the specific form of 
the stationary measures (5.4), it is clear that possible examples of monotone processes 
with stationary product measures of this form cannot be of misanthrope type. 


Lemma 5.1. A misanthrope process defined by the generator (5.1), that has stationary 
product measures and exhibits condensation is not monotone. 


Proof (5.5) gives necessary conditions for the monotonicity of the misanthrope process 
and implies with (5.4) that 

w{n - 1) ^ r(n,0) 

w{n) r(l, n — 1) 

is non-decreasing. This implies that the ratio converges to fc £ (0, oo], which is the reg¬ 
ularity assumption (2.10). Assuming the process condenses in the sense of Definition 
2.1, then Proposition 2.3 implies fc < oo. Now we have 

^ 4>c win) ^ w{n - (5.7) 

w[n) 

for all n G N. Therefore, w{n) ^ w{0)(pf^ which implies 

N N 

w{n)4>'2 ^ rt;(0) —)• oo as —)• oo . (5.8) 

n=0 n=0 

We conclude that the critical partition function diverges and the critical measure 
does not exist, which is a necessary condition for condensation. Therefore condensation 
does not occur in misanthrope processes with stationary product measures. □ 
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In [17] generalised misanthrope processes have been introduced where more than 
one particle is allowed to jump simultaneously. They are defined via transitions rj —)• 
r]+n{5y—5x) forn G {0,..., 'qx} at rate {y—x) and conditions on the jump rates for 
monotonicity are characterized. This class provides candidates for possible monotone, 
condensing processes with product measures as we discuss in the next subsection. 

5.2. Generalised zero-range processes. The generalised zero-range process (gZRP) 
[17] is a stochastic particle system on the state space defined by the generator 

C^Z^Pf{r^)= '^ak{Vx)p{x,y){f{v''^^’''>^) - f{q)) . (5.9) 

x^yGA k=l 


Here € Ql is the configuration after k particles have jumped from x to y G A. 

The jump rates ak{n) satisfy ak{n) = 0 if k > n, and we use the convention that 
empty summations are zero. We consider translation invariant p{x, y) on a finite lattice 
A = {1,..., L} and note that the process preserves particle number Yhx Vx = 

It is known [27,43] that these processes exhibit stationary product measures if and 
only if the jump rates have the explicit form 

akin) = , ( 5 . 10 ) 

n[n) 

where ( 7 , /i : N —)• [0, 00 ) are arbitrary non-negative functions with h strictly positive. 
The stationary weights are then given by w{n) = h{n). Monotonicity of the gZRP can 
be characterized in terms of 


n—k 

^ ^ (*^n—m(^) Q^n+1—“t“ 1)) ■ (5.11) 

m=0 

The gZRP is monotone if and only if 

Rkin) ^ 0 for all n ^ 1 and A: G {1,..., n} 
akin + 1) ^ Rkin) for all n ^ 1 and /c G {1,..., n} . (5.12) 

We note these conditions arise from a special case of the results in [17, Theorem 2.11] 
on generalised misanthrope models, since akin) depends only on the occupation of the 
exit site and not the entry site. 

In this class, which is also discussed in detail in [27], condensing processes which 
are monotone, homogeneous, and have stationary product measures with a power tail 
win) ~ n~^ with b G (1,3/2] are conjectured to exist. As an example, consider the 
gZRP with rates given by 


akin) 


'o 

< k-\f - ^)-^ 


if fe = 0 or n = 0 
if A: G {1,..., n — 1} 
otherwise . 


(5.13) 




MONOTONICITY AND CONDENSATION IN HOMOGENEOUS STOCHASTIC PARTICLE SYSTEMS 


25 


Since ak{n) is of the form (5.10) the process exhibits stationary product measures with 
weights of the form 


1 if n = 0 
n~^ otherwise 


For all 6 > 1 and L ^ 2 the ratio converges to ^ as ^ oo [9] so by 

Proposition 2.3 the process exhibits condensation. To prove the process is monotone 
we must show the rates satisfy the conditions given in equation (5.12). We first prove 
Rk{n) ^ 0 for all S {1,..., n} and n ^ 1. Since a„(n) — Q;„+i(n) = 0 for all n ^ 1 we 
can drop the m = 0 term from the definition of (n). We have 


n—k 


Rk{n) = 


m 


m=l 


^ m 
n 


-b 


- 1 - 


m 


n + 1 


-b-i 


Since (1 — x) ^ is increasing in a; for 6 > 0 and m/n> m/{n + l) we have 


Rk{n) > 0 for all /c € {1,..., n} and n ^ 1 . 


We also need to show ak{n + 1) ^ Rkin) for all k e {1,... ,n} and n ^ 1. Taking 
discrete derivatives in k 


ak+i{n + 1) - Rk+i{n) - {ak{n + 1) - Rk{n)) = ak{n) - ak{n + 1) 

-b 


= k- 


1 - 


k 


n 


-k- 


1 - 


k 


n + 1 


> 0 , 


so ak{n + 1) — Rk{n) is an increasing function in k. Therefore, 

Oikin + 1) — Rk{n) ^ Oii{n + 1) — iii(n) 
for all k ^ 1, and it suffices to show 

A{n) := Q!i(n + 1) — Ri{n) ^ 0 for all n ^ 1 . (5.14) 

We present numerical evidence in Figure 3 which corroborates our claim that the pro¬ 
cess with rates (5.13) is indeed monotone for b € (1,3/2] and is not for b > 3/2. 


5.3. Homogeneous monotone processes without product measures. The chipping 
model is a stochastic particle system on the state space Ql = introduced in [44,45]. 
The dynamics are defined by the generator 

Y >0)p(x,y)(/(7?"-^)-/(r?)) 

x,y£AL 

+ Y > 0)p{x,y){f{v + r]xiSy - ^x)) - fiv)) ■ (5.15) 

x,y&AL 

Here ij + rix{6y — 5x) denotes the configuration after all the particles at site x have 
jumped to site y, which occurs at rate 1, and single particles jump at rate re > 0. The 
spatial part p{x, y) is again spatially homogeneous as described in Section 5.1. 

It is easy to see that a basic coupling will preserve the partial order on the state space 
VLl as defined in Section 2.2. Therefore, by Strassen’s theorem [39], the chipping model 
is a monotone process and Lemma 2.2 implies that conditional stationary measures of 
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n 


Figure 3. Monotonicity condition (5.14) for b = 1.25, b = 1.5 and 
b = 1.65. For b = 1.65 the function A{n) falls below zero, implying the 
gZRP with rates (5.13) is non-monotone. For b = 1.25 and b = 1.5 the 
function A{n) is positive indicating the process is monotone. 


the process are ordered in N. The condensation transition in the chipping model was 
established on a heuristic level in [44-46]. We have defined the critical density pc 
only for systems with product stationary measures (see (2.15)). In general, the critical 
density on a fixed system of size L ^ 2, with unique invariant measures pl,n, can be 
defined as the background density of bulk sites 

Pc{L) := limsup ^ (5.16) 

N^OO ^ ^ 


Notice if are conditional product measures (see (2.9)) then pc{L) is consistent 
with (2.15) and in particular independent of L, which follows from Proposition 2.3 
(more explicitly (2.12)). For the chipping model in the case L = 2, the process re¬ 
duces to a 1-dimensional process on {0,..., A^} and the measure p 2 ,N and pc{2) can be 
computed explicitly to find 

P^{2) = ^(V2^FTT-1) . (5.17) 

In [44-46] the critical density in the thermodynamic limit is defined as 

Pc := sup |/9 ^ 0 : —)• 0 as L —)• oo such that ’ 

inspired by the fact that in case of condensation the second moment is either dominated 
by the condensate and scales like the system size L, or it diverges since the maximal 
invariant measure does not have finite second moment. It is shown by heuristic com¬ 
putations in a mean-field limit that 

Pc = Vw + 1 - 1 . 


This suggests that the critical density can depend on the system size L for distributions 
with non-product stationary measures. 
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The yTe scaling of the critical density can be intuitively understood in the two site 
chipping model with N particles. This process can be interpreted as a symmetric ran¬ 
dom walk on the state space {0,..., with jumps i —)• i ± 1 at rate w and random 
jumps to either boundary (resetting, i —^ 0 or A^) at rate 1. After a reset the particle 
diffuses at rate w and reaches a typical distance yTc from the boundary until the next 
reset. So this model is a monotone and spatially homogeneous process that heuris- 
tically exhibits a condensation transition with finite (size dependent) critical density, 
but it does not exhibit stationary product measures. Condensation is also observed in 
models where chipping is absent (w = 0) and the dynamics result in a single block of 
particles jumping on the lattice {1,..., L} corresponding to the critical density pc = 0. 


Appendix A. Connection to statistical mechanics 


Condensation and non-monotonicity are also related to convexity properties of the 
entropy, which we briefly describe in the following in a non self-contained and non¬ 
rigorous discussion that is aimed at readers with a background in statistical mechanics. 
In the thermodynamic limit the canonical entropy is defined as 

s(p) := lim y log Zl^jv . (A.l) 

Z/—^CX) / V 
N/L^p 

For the processes we consider, equivalence of canonical and grand-canonical ensembles 
has been established in [10] for condensing or non-condensing systems, so s{p) is given 
by the (logarithmic) Legendre transform of the pressure 


p{(j)) := logz((/>) . 


(A.2) 


This takes a particularly simple form since the grand-canonical measures are factoris- 
able, and is a strictly convex function for 4> < (pc- General results then imply that s{p) 
also has to be strictly convex below the critical density pc (see e.g. [47]), which holds 
for non-condensing systems and condensing systems with pc = oo. For condensing 
systems with finite critical density s{p) is linear for p > pc, consistent with phase sep¬ 
aration phenomena, where in this case the condensed phase formally exhibits density 
oo (see e.g. [32] for a general discussion). 

It is not possible to derive general results for finite L and N, but if we assume that 
the ratio of weights w{n — l)/w{n) is monotone increasing in n, we can show that a 
monotone order of 7rL,N implies that N logZ^^M is necessarily convex. Note that 
with (2.10) our assumption implies that w{n) has exponential tails with (^c G (0, oo) 
or decays super-exponentially with (f)c = oo, and in both cases the system does not 
exhibit condensation. We can define rcj—1) = 0 so that w{rix — ^)/'w{'>lx) is a monotone 
increasing test function on (1^. It is easy to see that for its canonical expectation we 
have for all L ^ 2 and N ^ 2 


T(L,N 


wjrjx - 1 ) \ 
w{rix) ( 


Zl,n-i 

Zl,n 


(A.3) 


Therefore, monotonicity of the canonical measures implies that the ratio of partition 
functions (A.3) is increasing and the discrete derivative of log Zl^m in N is decreasing. 
We expect that in the limit L —^ oo the monotonicity assumption on w{n — 1)/w{n) is 
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not necessary, and ^ log Z/, at is convex in N for all non-condensing systems, consistent 
with Strict convexity of s{p). 

For condensing systems the weights w decay sub-exponentially, and if w{n — l)/w{n) 
is monotone then it has to be decreasing in n. Therefore the choice u;(—1) = 0 im¬ 
plies /(ry) = w{rjx — l)/w{rjx) is not a monotone test function, and the above general 
arguments cannot be used to relate non-convexity of log to the absence of a 
monotone order in For particular condensing systems, however, it has been 

shown that ^ log is typically convex for small N < pcL and concave for larger 
N > PcL [12,25]. These results focus on power law and stretched exponential tails 
for w{n), and have been derived for zero-range processes where the ratio tv 

is equal to the canonical current. Non-monotone behaviour around the critical density 
therefore has implications for finite-size corrections and derivations of hydrodynamic 
limits as mentioned in the introduction. 


Appendix B. The infinite mean power law case with b = 2 


Consider stationary weights of the form w{n) = n~‘^ with w{0) = 1, we prove the 
non-monotonicity of Zl,n/'w{N) in a similar fashion to the proof of Lemma 3.6 sum¬ 
marised in the following lemma. 


Lemma B.l. For stationary weights of the form w{n) = n ^ with w{0) = 1 we have 


N 


Z2,N 
log{N) \w{N) 


— 2z{l) ) F 2 = 4: as N ^ oo . 


For L > 2we have 


N 


Zl. 


N 


- Lz{l) 


L-l 


Fl := z{1)Fl-i + {L - l)z{l)^ ^^2 as N ^ 00 , 


\og{N) \w{N) 
which is positive for all L ^ 2 since F 2 > 0. 


(B.l) 


Proof First consider the case L = 2. As in the proof of Lemma 3.6 we will utilise the 
full Taylor expansion of (1 — x) integral bounds on monotone series, and assume N 
is even, for N odd there exists obvious modifications to the proof. We have 


Z2,N 

w{N) 


N /2 

2z{l) = 

n=l 



n=l 


- 2‘^N-‘^ . 


Where the terms n = 0 in the above summations cancel. Substituting the Taylor expan¬ 
sion of (1 — we find 


y 00 N/ 2 cx:) 

^^^-2^(1) = 2^(i + l)Ar-*^n"2+*-2 ^ n“2-2^iV-2. (B.2) 

i=l n=l n = N / 2 +l 


Now we are in a position to apply the integral bounds (3.16), first noting that 
is decreasing for i = 1, constant and equal to 1 for i = 2, and increasing for i ^ 3. 
Multiplying both sides of (B.2) and applying the integral bounds it is easy to show 


N / Z2,n 
log(A^) \w{N) 



—)• 4 as A" —)• 00 . 
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Now consider the case L > 2 and make the following inductive hypothesis 


lim 


N [ Zl,n 


- Lz{l)^-^ ]=Fl = zil)FL-i + (L - 1)z(1)^-"F2 


L-2 i 


IV-J-oo log(A^) \w{N) 

As in the proof of Lemma 3.6 write 

N ( -^L+l.jV 
log(A^) \ w{N') 


-(i +1)2(1)' 


^ (’F'^ ui(lV-n) 

^ n=0 ^ ' 


log(iV) 


-^(1)^ + 


N 


. V/2-1 


log(iV) V ^0 


( Y1 


ZL,N-n 

w{N) 


-Lziiy 


^L,N 


'Jl,n 


(B.3) 


We fist establish the limit of Ql,n in (B.3). The inductive hypothesis can be rewritten 
as 




w[n] 


= [Fl + on 


log(iV) 

N 


+ Lz(l) 


L-l 


(B.4) 


Similar to the proof of Lemma 3.6 Ql,n can be written in the form 

^V/2-l 


Ql,n = 


N 


log(iV) 


Fl + ON 


+ Lz{l) 


L-l 


N 


log(iV) 


w{n 

n=0 
N/2-1 
yy w{n 

n=0 


w{N — n) log(A^ — n) 
w{N) N — n 

w{N — n) 


w{N) 


-z{l) 


(B.5) 


Since log(A^ — n) is decreasing for n G {0,... N/2 — 1} we can find upper and lower 
bounds of the first term, by pulling out the logarithm, of the form 


log(iV/2 - 1) 
log(iV) 

N 

log(iV) 


Fl + ON 


'N/2-1 


w{n) 




Fl + ON 


n=0 

''n/2-1 

Y 

n=0 
^ N/2-1 


^ w{N — n) N 
w{N) N — n 


^ (.Fl + ov(1)) Y^ 


w{N — n) log(A^ — n) 
w{N) N — n 

w{N — n) N 


n=0 


w{N) N — n 


Applying the same steps exactly as they appear in the proof of Lemma 3.6 to the upper 
and lower bounds we have 
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To identify the limit of at in (B.3) we again follow the steps given in the proof of 
Lemma 3.6, which implies 

lim El,n = lLzil)^-^F2 . (B.7) 

A'—>-oo Z 

Combining this with (B. 6 ) we have 

^ as iV ^ oo . 

From the recursion (B.l) it is obvious that Fi will have the same sign as F 2 , completing 
the proof of Lemma B.l. 

□ 


Appendix C. On the sign of ^2(6) 

In this section, we compute the sign of ^ 2 ( 6 ) for 6 G (1,2), where 


00 ^ 2 — 1 


= 2^:1 no +'’>1 


'ib—l—i 


^b-1 


-2- 


*=i-i=o --^ + * 

Recall the definition of the Pochhammer symbol 

1 if n = 0 

{q){q + 1)... {q + n — 2){q + n — 1) for n ^ 1 

and the hypergeometric function 


iQ)n = 


2 Fi{c,d,e,z) = V - 

^ ^ 7 


F- {c)i{d)i 


i=0 


(e)i 


We now show 


^ 2 (^) = -^ 2 Fi (l-b,b,2-b,^ 


(C.l) 


which in particular implies ^ 2 ( 8 / 2 ) = 0 by evaluating the hypergeometric formula. 
Factorising the term 2^/{b — 1) from ^ 2 ( 6 ) and rearranging terms inside the summation 
we have 


00 ^ / X ^ 




2 = 1 


2 ^ A 1 /I 


5-1 


5- 


1=0 
2—1 




— b + i 
5-1 


- 1 


1 — 5 + i 


i=0 ^ 2 j—Q 

Now use the following identities to simplify the terms inside the summation 


2—1 


(2 - 5), 


YlU + b) = {b)i and (1 - 5 + i) = (1 - 5) 


which gives the required result (C.l). 
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To complete the proof we use the following two relations for hypergeometric func¬ 
tions, Euler’s transform [49, 15.3.3] 


2 Fi{c,d,e,z) = (1 - zY '' 2-^1 (e - c,e - d,e,z) , 

and Gauss’s second summation theorem [49, 15.1.24] 


e—d—c 



Therefore, 



0 F 22 '>-ir (2 - b) 


To calculate the sign of ^ 2 ( 6 ) we first note that the gamma function r( 3 :) is positive 
for all a; > 0 and negative in the region — 1 < a; < 0 , which implies 



< 0 for 6 e (1,3/2) 
> 0 for 6 e (3/2,2) 
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